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ABSTRACT 

Temperature data, collected over the past 40 years, was 
compiled and averaged in the first known study of mean ther- 
mal conditions throughout the waters of Monterey Bay. The 
following results were obtained: 

(1) The distribution of mean sea surface temperature 
in the Bay was obtained by calculating mean monthly values 
at selected grid points and drawing isotherms, 

(2) Progressive warming of the upper 100m was observed 
to occur in a nearly linear fashion with time, resulting in 
a maximum temperature increase of 1.8°C over the 29 year 
interval from 1931 to 1960. This implies that the in- 
tensity of upwelling is also decreasing v.'ith time, 

(3) Geographic variations in the rates of upwelling 
and downwelling, causing a relatively "warm tongue" along 
the Canyon axis, appeared to be a major element influencing 
the dynamics of the Bay, and 

(4) Mean currents were inferred from the topographies 
of mean monthly isothermal surfaces. Limited comparisons of 
the inferred flow to past measured currents were very en- 
couraging and tended to support the feasibility of such an 
inference in the shallow waters and irregular topography of 
Monterey Bay. 
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I. INTRODUCTION 



Monterey Bay is a large, semi-elliptical indentation 
in the California coastline separated into nearly equal 
halves by a deep submarine canyon (Figure 1) . ^ The outer 
extremity of the Bay may be defined by a line joining Point 
Santa Cruz, to the north and Point Pinos to the south. 

A. THERMAL CONDITIONS IN MONTEREY BAY 

Several studies have been conducted of thermal condi- 
tions at specific locations in Monterey Bay. The most ex- 
tensive were those by Skogsberg [1936] and Bolin [1964] , 
both spanning five year periods. Skogsberg studied the 
southern portion of the Bay from 1929-1933 and Bolin re- 
searched conditions at a station in the Monterey Submarine 
Canyon (see Appendix A) . Although these studies \vell de- 
fined the thermal conditions for a portion of the Bay, ap- 
parently no one had investigated the Bay as a whole. 

As discussed in Appendix A, the author com.piled a 
considerable amount of temperature (and some salinity) data 
covering the period of time from Skogsberg 's study until 
the present. As outlined by the author's objectives (Para- 
graph C) , this data was utilized to study mean monthly ther- 
mal conditions throughout the Bay and to investigate the 
possible relationship between the thermal conditions and 
indicated m.ean currents. 

1 

All figures will be found in Appendix C. 
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B. 



CURRENTS IN MONTEREY BAY 



1 . Importance 

The distribution of current in Monterey Bay is 
virtually unknown but there is little doubt of its great 
importance in determining the ecological future of the area. 
Pollutants, transported by the currents, are scattered 
freely along the beaches in the form of sewage and general 
rubbish. Less evident is the adverse effect on marine life 
created by the discharge of industrial pollutants into the 
Bay. Currents not only distribute "harmful'' pollutants but 
also serve to maintain the chemical and thermal balance 
necessary to support marine life. 

2 . Why Little is Known Concerning Distribution of Bay 

Currents 

The well established currents of the open oceans 
are relatively strong and unchanging with time (e.g., the 
Gulf Strear.i and the offshore California Current) and hence, 
their general movements are fairly well understood. Monterey 
Bay, being enclosed on three sides, presents a different 
situation. The Bay currents are weak and variable and thei r 
mean flow is almost completely disguised by the effects of 
local winds, tides, topography, and seasonal changes in the 
offshore currents. 

Past current studies in the Bay have been of a 
limited nature, both in time and in geographical coverage. 
They have contributed towards an understanding of "short 
term" currents (such as those resulting from the tides) but 
too few observations have been made to be of much use in 
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describing the mean flow. In time, and with many addi- 
tional current observations, a study of the combined ef- 
fects could prove to be of great value. 

3 . The Feasibility of Approximating Mean Currents by 
the Geostrophic Approach 

Indirect approximations of current using the geo- 
strophic theory may be divided into two categories: (1) 

the classical use of the geostrophic method, and (2) ap- 
proximations derived from the classical method. 

a. The Classical Geostrophic Method 

Basically, flow obtained by classical applica- 
tion of the geostrophic method is that which results from 
the balancing of the Coriolis force and the pressure force 
between ti^o points on a surface of variable elevation, 
(Neumann and Pierson [1966] discuss the theory in their 
book entitled. Principles of Physical Oceanography ) . Com- 
puted elevation of the sea surface at a given location in 
the ocean is represented by a term called the "dynamic 
height," a relative parameter based on the density of sea 
water at the location and referred to an assumed level-of- 
no-motion. A high value of density (usually colder water/ 
yields a lower dynamic height and conversely. By geo- 
strophic theory, the circulation in the Northern Hemisphere 
is counterclockwise about locations of lower dynamic height 
(depressed surfaces) and clockwise about locations of higher 
dynamic heights (elevated surfaces) . The level-of -no-motion 
is the depth of water (as best can be ascertained) where 
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absolute current velocity equals zero and this rarely oc- 
curs at depths less than 150m. 

Accurate representation of the mean current 
based on the classical geostrophic method thus requires two 
prerequisites; (1) a knowledge of mean temperature and 
salinity (to obtain mean density) o£ the water column at 
several well spaced locations in the area of interest, and 
(2) a level-of -no-motion everywhere shallower than the 
depth of the bottom (otherwise, there is no reference ex- 
cept the bottom itself from which to compute true velocity) . 
Finally, of course, wind driven and tidal currents are 
neglected in computing means. 

b. Approximations Derived from the Classical 

Method 

Several approximations to the classical ap- 
proach have been used in relatively deep waters with con- 
siderable success. Barker and Denham [1970] used the 
classical method to plot geopotential topographies (sur- 
faces of equal dynamic height) for various levels-of-no- 
motion off the northeast coast of New Zealand. The close 
correlation between the shapes of the topographies sug- 
gested that broad current patterns could be approximated 
without using an accurate level- of -no-motion . Duncan and 
Nell [1969] demonstrated that isopycnal (equal density) 
surfaces were closely related to current patterns in their 
study off the Cape Coast. Drift cards were employed to 
measure general current trends and were found to follow the 
isopycnal contours as predicted by geostrophic theory for 
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the Southern Hemisphere. Leipper [1970] demonstrated that 
temperature data alone could be used to approximate cur- 
rent in the Gulf of Mexico. Monthly topographies of the 
22°C isothermal surface were constructed by him from data 
collected on various cruises during the period 1965-1966. 

The current was found to closely follow the contours, moving 
in a direction placing colder waters to the left of the 
flow. These results were verified by geostrophic flows com- 
puted at various locations in the Gulf. This approach is 
based on the assumption that temperature changes are di- 
rectly proportional to density changes throughout the vi^ater 
column, and thus isothermal surfaces may be considered in 
place of isopycnal surfaces. 

Application of any of the above approximations 
to the shallow waters and irregular topography of Monterey 
Bay would be a considerable extension of usual practice. 
However, the theory of geostrophic flow does not depend on 
the depth of water, therefore, geostrophic methods might be 
useful in indicating broad current patterns in the shallov; 
waters of the Bay. Considerable temperature data was com- 
piled for Monterey Bay in this study (see Appendix A) . The 
author used this data to investigate the feasibility of 
relating mean currents to mean isothermal surfaces in a 
manner similar to that of Leipper. 

C. OBJECTIVES 

The objectives of this study are to: 
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1. Compile temperature and salinity data from previous 
oceanographic studies conducted in Monterey Bay, 

2. Conduct a study of the distribution of mean monthly 
temperature in Monterey Bay and construct charts of the mean 
monthly distribution of sea surface temperature, 

3. Investigate the feasibility of inferring mean cur- 
rents from the topographies of mean isothermal surfaces, and 

4. Prepare a computer program for the continuing 
storage and retrieval of oceanographic data collected in 
Monterey Bay. 

D. ORGANIZING AND USING THE DATA 

The IBM 360 computer at the Postgraduate School was 
used to compile the temperature data and to compute monthly 
averages (as discussed in Appendix B) . A grid system, 
dividing the Bay into 19 blocks, 4mi by 4mi in dimensions, 
was used to segregate the temperature data (Figure 1) . The 
author's initial grid system covered a larger area (Figure 2) 
containing 39 blocks, but only the numbered blocks held 
significant data. The mean temperature of a given block, 
based on the averaging of all temperatures compiled for that 
block, was assumed to be effectively located at its center. 

The temperature data was compiled over the 40 year 
interval from 1929 through 1968. Temperatures collected by 
all reliable methods were combined in the averaging process 
(although the computer program did allow for data to be 
retrieved by type, i.e., Nansen, XBT, BT) . 
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The risks inherent in the "lumping” together o£ all 
types of temperature data spanning a 40 year interval and 
further reducing data from 4mi by 4ni blocks to represent 
points in their centers were fully realized. This lumping 
was essential due to the erratic distribution of data within 
the Bay and also because no reasonably short time span held 
enough data to support this study. The errors introduced 
by the above procedures are discussed in Chapter II. 
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II. VALIDITY OF THE COMPUTED MEAN MONTHLY TEMPERATURES 
(AND DENSITIES) 



A discussion o£ the questions raised in the Introduc- 
tion involving the method of computing mean monthly temper- 
atures is considered necessary before an appraisal of the 
results (Chapters III and IV) can be made. 



A. TIME RANGE OF DATA 

Data from three separate time intervals within the 
previously mentioned 40 year time span was averaged and 
compared to determine the degree of consistency of thermal 
conditions relative to the progression of time. Three Bay 
studies (as referenced in Appendix A) were used: (1) 

_ f -I r\ r\ /'irvri ^r»rr^ 

i. ^ j o j , xjOlxii X o X “ X ^ ^ j 

CalCOFI surveys (1954-1967).^ The time intervals between 
midpoints (mean years) of the studies were 22 years from 
Skogsberg to Bolin and 7 years from Bolin to CalCOFI , with 
a total time interval of 29 years (1931 to 1960) . 

The annual thermal cycles of the three periods (Fig- 
ures 4-6) were generally similar in shape with upwelling 
noticeable from the surface to 500m (as characterized by 
the humps in the isotherms during the months of April 
through June) . A consistent feature over the 29 year in- 
terval was the apparent weakening in intensity of upwelling 



2 

The oceanographic stations used in these studies are 
shown in Figure 3. Skogsberg 's station is denoted by "C’’ , 
Bolin's by "B" and CalCOFI by a square numbered "5". It is 
noted that the three stations are contained within a radius 
of tv.'o mi. 
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with time. The 9°C isotherms from all three studies are 
plotted in Figure 4 to demonstrate the effect. It was ob- 
served that the depths of the 9°C isotherms were relatively 
close together during the winter months. During the re- 
mainder of the year (especially noticeable in May) the 9°C 
isotherms were observed to deviate from one another. The 
deviation followed a consistent pattern with the 9°C iso- 
therm of the earliest study (Skogsberg in 1931) at the 
shallowest depth and the 9°C isotherms of the later two 
studies at greater depths as time progressed. This "flat- 
tening out" of the 9°C isotherm with time could indicate a 
weakening in intensity of upwelling. 

A general warming with time was also observed to occur 
throughout the waters of the Bay. Figure 6A depicts the 
temperature increase (AT) observed to occur over the period 
from 1931 (Skogsberg' s mean) to 1960 (CalCOFI's mean) at 
selected depths throughout the year. The curves represent 
the variation of AT for each month and for depths from the 
surface to 500m. Figure 6A shows that the maximum tempera- 
ture increase occurred at depths from the surface to 100m 
and during the period from October to December. At depths 
from 200-500m, AT cycled between positive and negative 
values throughout the year; however, AT was always positive 
(temperature increase) in the months of May, June, September, 
and October. 

The consistently positive values of AT in May and June 
for depths from the surface to 500m might be another indica- 
tion that upwelling was decreasing in intensity (with time) 
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throughout this region, since the warming in May and June 
would also indicate a flattening out of the isotherms. 

The reasons for the maximum temperature increase during the 
months of October to December are unknown. It could be re- 
lated to changes in the "Oceanic Period," or possibly, to 
changes in the transition period between the influence of 
the California and Davidson current systems. 

In summary, the general warming of the water column 
with time was observed to be consistent in the upper 100m 
(by the similarity in the shape of the AT curves in Figure 
6A at these depths). At depths below 200m, the consistency 
was not as good (as observed by the AT curves cycling be- 
tween positive and negative values); however, the total 
range of AT was only 0.5°C in the deeper waters. 

The maximum temperature increase over the 29 year 
interval was 1.8°C (the 50m curve in November from Figure 
6A) resulting in about l/15th of a degree C increase per 
year. This small annual temperature change, along with 
the general similarity of the curves (as observed in Fig- 
ures 4-6) was considered to represent consistent behaviour 
with time. It was concluded that the author's method of 
averaging temperature values was justified. This provided 
averages that were representative of the overall period. 

B. POSITIONING OF STATIONS 

The computer grid system was based on the assumption 
that all temperature values found in a given block were lo- 
cated at its center. The grid system was initially positioned 
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so as to minimize the distance of known oceanographic sta- 
tions from grid centers. The maximum deviation from such 
centers was 1.6 mi. Data was assumed to be scattered 
randomly around the centers of the blocks in the grid sys- 
tem. 



C. USE OF THREE DIFFERENT DATA COLLECTION SYSTEMS 

The data was collected by Nansen, BT , and XBT devices. 
Several papers have been written concerning system errors 
peculiar to each. The magnitude of the errors in tlie BT ' s 
and XBT ' s are apparently undeterminable. 

To investigate the significance of possible system 
errors, the standard deviations (a) of temperature (and 
density) values were included in the computer output. 

(X ± a) represented the interval containing 67% of all 
values averaged for a given mean (X) . 

Standard deviations of mean monthly sea surface temper- 
ature versus time (by months) were plotted for selected 
blocks. Figure 7 (Block 13) represented the worst condi- 
tion found. The larger standard deviation in July was not 
representative; the consistent features of all graphs were: 

(1) minimum values of standard deviation for the month of 
January, and (2) the relatively small variation of the 
values throughout the year. 

Plots of standard deviation versus depth (Figure 8) 
were also consistent. Values of standard deviation halved 
within 30-40m of the surface, reached values of 0.5°C anywhere 
from 50-275m and then linearly decreased to minimum values of 
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0.1-0.2°C in deeper waters. The decrease of standard devia- 
tion with increasing depth of water and its practically 
invariable value (0.1-0.2°C) below 800m, where the data was 
exclusively from Nansen casts, indicated that seasonal 
changes were mainly responsible for variations in mean 
temperatures. System errors were thus neglected. 

D. GEOGR_\PHIC DISTRIBUTION OF DATA 

The 3000 oceanographic stations used in this study 
(25,100 temperatures) were erratically distributed through- 
out the blocks in the computer grid system (Figure 2) . 

Blocks 1-6, 13, 20, and 21S contained the major portion of 
the data (10-85 surface temperatures per block per month) 
and thus were given major consideration in the analysis. 

The remainder of the blocks contained random amounts of 
data, except for Blocks 17 and 19 which generally contained 
3-6 values per month each. 

Clearly, the blocks containing relatively few tempera- 
tures (less than 10 per month) would yield less reliable 
averages than blocks with many values. Choosing a minimum 
number of temperature values to achieve relatively accurate 
mean monthly temperatures was based on a study of continuity 
between weak blocks (with less than 10 temperatures per 
month at a given depth) and strong blocks (more than 10 tem- 
peratures) . 

Graphs of annual mean temperature cycles, m.ean monthly 
temperature sections, and standard deviation versus time and 
depth ivere drawn for weak Blocks 8, 17, and 19 and for strong 
Blocks 1-6 and 13. The annual m.ean temperature cycles of the 
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weak blocks compared closely to those o£ the strong blocks. 
This comparison \-ias well represented by Block 19 (Figure 9) 
and Block 3 (Figure 10) . 

Further continuity in data between weak and strong 
blocks was demonstrated by Figure 11, showing the thermal 
structures of two vertical sections (see Figure 12 for lo- 
cation of sections) , one oriented north to south connecting 
Blocks Ij 2, 17, and 5 (S-1) and the other oriented along 
the Monterey Submarine Canyon axis (S-2) connecting Blocks 
13, 3, 17, and 19. The isotherms were smooth and continuous 
in the upper and lower waters for both sections throughout 
the year. In depths from 50-100m, the isotherms were smooth 
and continuous along section S-2. Along section S-1 the iso- 
therms were occasionally depressed at Block 17. The depres- 
sion of isotherms at Block 17 in section S-1 and lack of 
depression for the same block in section S-2 will be shown 
to be consistent in Chapter III, paragraph B. 

Finally, the standard deviations of mean monthly tem- 
peratures versus time and depth were comparable for all 
blocks, weak or strong (Figure 7 and 8 were representative). 

The continuity exhibited in the vertical sections be- 
tween adjacent weak and strong blocks (and within the blocks 
as time progressed) supported the limited use made of the 
weak blocks for constructing mean monthly isothermal sur- 
faces in Chapters III and IV. As the isothermal surfaces 
were prepared, the weak blocks demonstrated further consis- 
tency of data, especially Blocks 8, 17, and 19. Examples of 
this are noted in Chapter IV. 
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III. THERI'IAL CONDITIONS 



A. THE ANNUAL THERMAL CYCLE IN MONTEREY BAY 

The "3-phase” annual thermal cycle occurring at specific 
locations in Monterey Bay has been described by Skogsberg 
[1936], Bolin [1964] and others. The following discussion 
considers the variation in mean thermal conditions relative 
to various geographic locations throughout the Bay, a sub- 
ject apparently not previously dealt with. 

Figure 12 depicts three geographic regions in the Bay 
that the author will refer to in this Chapter. The choice 
of the 50 fathom curve as a dividing line was arbitrary, but 
gave the desired sectioning. 

1 T’'U o P o-nAT-or^ P o rr t on 

Since Bolin described the thermal conditions for 
the Canyon region a detailed account will not be repeated 
here; however, a brief description of the "phases" will be 
necessary for later comparisons. 

a. The "Upwelling Period" 

Commonly considered to span from mid-January to 
September, the Upwelling Period is characterized by the 
rising of the deep, colder waters. As found in this study 
(Figure 10) , the cold waters began their upward movement in 
mid- January from depths as great as 700m (however, the move- 
ments were minimal below 300m) . The maximum rate of rise 
occurred between February and March and the peak was found 
in April or May when the 11°C water nearly reached the 
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surface. Although upwelling reached its peak later at 
greater depths, none was observed after July. 

The descent (downwelling) of the cold waters 
generally commenced in June and was apparent until late 
November at greater depths. 

b. The "Oceanic Period" 

This period is commonly considered to occur 
during September and October and is i\:ell defined by the 
sharp vertical temperature gradients at depths from the 
surface to about 100m. 

c. The "Davidson Current Period" 

In November the south-moving California Cur- 
rent is displaced by the north-moving Davidson Current along 
the coast. The influx of warmer waters into the Bay re- 
sults in a weakening of the vertical temperature gradients 
and a resulting well-mixed layer of water from the surface 
to about 50m. The northerly flow continues until mid- 
January when the California Current again predominates and 
the "3-phase" annual cycle commences again. 

2 . The Regions of the Shallows (see Figure 12) 

The annual thermal cycle in the Shallows was simi- 
lar to that in the Canyon Region with two noteworthy ex- 
ceptions : 

First, the very shallow regions at the extreme 
northern and southern ends of the Bay possessed unique 
annual thermal cycles at depths from the surface to 10m. 
Referring to Figure 13, the annual cycle of the surface 
temperature in the Canyon (Block 3) showed a clear decrease 
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in temperature from January until early May and then a 
steady increase until late September (this was typical for 
all the deep water areas) , The annual cycles of surface 
temperature for the very shallow waters (Blocks 1 and 5) 
showed no similar decrease in late winter but rather a 
steady increase from mid-January until late September 
(Figures 13 and 14). This steady temperature increase, 
unique for the shallow water regions (less than 75m) , was 
probably due to greater heating of the water column from 
bottom reflection and the insolation of nearby land masses. 
This greater heating of the shallow areas is clearly shown 
in the distribution of mean monthly sea surface temperatures 
as will be discussed in paragraph E (Figures 19-30) . The 
relatively deep regions of the Shallows areas (Blocks 2 and 
4) exhibited a very slight amount of cooling from January 
until May, but not nearly as much as did the canyon Region 
(see Figures 15 and 16) . The second exception is discussed 
below. 

B. THE VARIATIONS IN INTENSITY OF UPWELLING 

In studying the various graphs constructed for the 
analysis of thermal conditions, the author came upon a 
feature which could be an important element in the dynamics 
of the Bay in general. Figure 17 is introduced to show the 
months of the maximum rates of upwelling and downwelling 
during the year. Inspection of the annual cycle in Block 3 
shows maximum rates of upwelling from February to March, and 



25 




I 



maximum rates of downwelling (or heating) from July to 
August (note the slopes of the 10°C and 11°C isotherms). 

Vertical temiperature section S-1 (see Figure 12) shows 
a strange dipping of the isotherms at Block 17 for all of 
the above mentioned months (Figure 18 is representative) . 
Block 17 was a weak block (less than 10 temperatures per 
month) and thus the dipping was suspected. To attempt to 
verify the dipping of isotherms in the Canyon, two addi- 
tional sections were constructed (see Figure 12 for loca- 
tions) , both running across the Canyon axis (as did section 
S-1). Section S-3 connected Blocks 2, 3, and 4, and sec- 
tion S-4 connected Blocks 2, 3, and 6 (emphasis is placed 
on the fact that sections S-3 and S-4 contained only strong 
blocks with greater than 20 temperatures per month each) . ^ 
The results are clearly illustrated in Figure J8 which sho\vs 
marked dipping of the isotherms across the Canyon axis in 
all sections. This results in a relatively v^^arm tongue 
along the Canyon axis for the months of maximum rates of up- 
welling and downwelling. 

The dipping of the isotherms across the Canyon axis 
could have been caused by the up-slope entry of the cold up- 
welling waters into the regions of the Shallows. During 
February and March, when upwelling rates are highest, the 
cold waters are ascending from depths down to 300m. As the 



3 

Table V lists the number of temperatures compiled for 
the various blocks of the computer grid system. 
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cold, dense water spreads to the Canyon walls, its hori- 
zontal movement is interrupted, and thus it is accelerated 
up the walls and into the shallower adjacent regions. 
Conversely, in July and August, when downwelling rates are 
highest, the cold waters are descending from the near- 
surface regions. Their rates of descent are accelerated 
into the Canyon since the deepening topography offers no 
resistance to their sinking. The transition period (from 
April to June) between upwelling and downwelling showed no 
dipping of isotherms across the Canyon. 

No dipping of isotherms was observed along the Canyon 
axis. Thi.s v;ould be as expected since the depth of water 
along the axis (section S-2) was everywhere deeper than 300m 
thus affording no bottom influence upon the cold upwelling 
waters . 

The annual mean temperature cycles of the relatively 
deep regions (Figure 16) and the very shallow regions (Fig- 
ure 14) of the Shallows further demonstrated the up-slope 
surge idea. During the months of highest up/downwelling 
rates. Blocks 2 and 4 (adjacent to the Canyon) exhibited 
somewhat stronger upwelling (greater slopes of the 10°C 
and 11°C isotherms) than did Blocks 1 and 5 which were at 
the extreme northern and southern ends of the Bay. It thus 
appeared that the cold up-surging waters might have lost 
some of their momentum by the time they reached the very 
shallow regions of Blocks 1 and 5. 
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C. COMPARISON OF THF ANNUAL THERM/vL CYCLE IN THE BAY 

TO THAT OF AN OFFSHORE LOCATION 

Block 20, located 22 mi west- southwest o£ the Bay 
(Figure 2) was compared to Block 3 in the Canyon. Figure 
17 demonstrated two broad differences between the thermal 
conditions of the two locations. 

1 . Mixing in the Upper Layers 

The offshore station exhibited generally warmer 
surface waters throughout the year and better mixing of the 
upper waters (0-50m) from January to April. At onset of 
the "Oceanic Period" in September, the offshore and Bay 
stations both exhibited maximum vertical temperature gradi- 
ents . 

At the close of the "Oceanic Period" in November, 
the Bay station exhibited somewhat better mixing than did 
the offshore station, maintaining its more uniform surface 
waters until January. The weakening of the vertical tem- 
perature gradients in the Bay (as described by Bolin) was 
assumed to be the result of the influx of warm waters 
brought into the Bay by the north-flowing Davidson Current 
(during the period of November through January) . The more 
uniform mixing of the Bay waters might have been due to a 
greater influence by the Davidson Current nearshore. 

2 . Upwelling 

Upwelling is considered to result from the northerly 
winds deflecting the surface waters to the right and away 
from the coast, allowing the deeper cold waters to rise. The 
intensity of upwelling should then decrease as the distance 
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to the coast increases. Figure 17 demonstrated that this 
was indeed the case, Block 20 showing less intense upwelling 
than Block 3. 

D. PROGRESSIVE WARMING 

Chapter II, paragraph A, contains a discussion of the 
progressive warming found for the entire water column in 
the Bay over the 29 year interval from 1931 to 1960. It 
should be noted that Bolin [1964] found an opposite indica- 
tion of progressive cooling from 1951-1955. He investigated 
the ’’previous and subsequent” sea surface temperatures taken 
at selected locations along the west coast of North America 
and found that the years of his study were all relatively 
cool and that 1955 had the minimum temperature in a ’’long- 
term trend.” In Chapter II, evidence was presented indica- 
ting continual progressive warming over the 29 year interval 
noted above. It thus appears that rhe cooling observed by 
Bolin was a short-term cycle. 

E. THE DISTRIBUTION OF MEAN MONTHLY SEA SURFACE TEMPERATURE 

IN MONTEREY BAY 

Figures 19-30 depict the distribution of mean monthly 
sea surface temperature in Monterey Bay. 

1 . Description of Figures 

The dots on each figure (representing the centers 
of their respective blocks) served as the points from which 
isotherms were drawn. The values of mean monthly tempera- 
ture are directly below the points. Table I lists the total 
number of temperature measurements per block per month (at 
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TABLE I . The Range in the Total Number of Temperatures 
Averaged to Compute the Means for Blocks in 
the Bay 
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Note 1: Blocks 9-11, 14, and 18 seldom contained data. 

Note 2: Block 20 (see Figure 3) is located 22 mi offshore. 
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the surface and 50m) that were used to compute the mean 
temperatures. Mean temperatures are not listed (or used) 
in the figures if less than four values were available 
for averaging. Dashed lines represent probable extensions 
of the isotherms. 

The circle in the upper right hand corner of each 
figure displays the local mean wind condition (taken from 
the U. S. Naval Oceanographic Office "Pilot Charts," 
series H.O. 1401) and the mean offshore currents (taken 
from the U. S. Naval Oceanographic Office Publication H.O. 
570) . The local mean winds are shown as a series of ar- 
rows (fine lines) that fly with the winds. The lengths of 
the arrows represent the percentages of the month that wind 
were from the directions shown (the percentage shov>?n for 
the longest arrow serves as a reference) . The number of 
tails on a given arrow defines the wind force by the 
Beaufort Scale, as sho\>?n below; 



Number of Tails Mean Wind Force (kt) 
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The mean offshore current is displayed by an arrow (bold 
line) extending the diameter of the circle. Mean current 
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direction is with the arrow and its speed (miles per day) 
is displayed inside the inner circle. 

The upper left hand corner of each figure shows the 
mean sea surface temperature at a location 15 mi offshore. 
This data was taken from Naval Oceanographic Office Special 
Publication, SP-123, and represented mean monthly values 
based on ship's injection temperatures collected over the 
past 100 years. 

2 . Discussion 

a. The Annual Cycle of Mean Sea Surface Temperature 

The following two factors appear to influence 
the distribution of mean sea surface temperature. 

(1) Upwelling 

The coldest waters in the Bay occur during 
the period from February to May (with April being the coldest 
month of the year) coinciding with the upwelling period. 
Figures 20-22 show the intrusion of the colder waters from 
seaward. It appears that the main intrusion broadly oc- 
curs along the axis of the Submarine Canyon (note the 11°C 
and 12°C isotherms in Figure 22) . 

(2) The Variation in Mean Air Temperature 
According to data contained in the U. S. 

Naval Oceanographic Office Pilot Chart series, the mean air 
temperature is at a minimum in January, rises to an annual 
maximum in September to October, and then falls to a minimum 
again in January. 

The relation of mean sea surface to mean 
air temperatures is demonstrated by following the monthly 
position of the 13°C isotherm. In January and February its 
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absence implies generally colder waters. In March it was 
seen to intrude into the Bay near the Moss Landing area. 

In April, the 13°C isotherm appears to be pushed aside by 
the intrusion of the colder waters (Figure 22) . During the 
period from May to October, the 13°C isotherm was observed 
to move steadily seaward as the waters in the Bay continued 
to warm. Although the positions of the 13°C isotherms are 
nearly identical in September and October, the overall 
thermal distribution shows September to be the warmest month 
of the year (Figure 27). From October until January, the 
13°C isotherm was observed to move steadily towards shore 
indicating the broad cooling of the surface waters. 

b. Horizontal Sea Surface Temperature Gradients 
Horizontal temperature gradients are weakest 
during the period from November through January with nearly 
uniform temperature throughout the Bay (less than 0.5°C 
change) occurring in January. The uniform temperatures 
characterizing this period are probably due to the influx 
into the Bay of \\rarmer waters from the north-flowing 
Davidson Current (prevailing during the above period) . In 
February, the temperature gradients are stronger and con- 
tinue to increase in strength until reaching their annual 
maximum in July. The mean offshore sea surface temperature 
of 13.9°C in July indicates that strong gradients are also 
set up offshore as well as in the Bay. Steady weakening of 
temperature gradients occur from August until January, com- 
pleting the annual cycle. 
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It will be recalled that maximum vertical tem- 



perature gradients occur in the surface layers during the 
Oceanic Period from September to October (see Figure 10). 
Bolin related the occurrence of maximum vertical gradients 
to the relatively windless (and also the minimum velocity 
of the Cali.fornia Current) conditions prevalent during this 
period. The occurrence of the maximum horizontal sea sur- 
face temperature gradients in July (2-3 months earlier than 
the maximum vertical gradients) miglit be related to the 
greater velocity of the south- flowing California Current. 

The cooler waters carried into the Bay by this current move 
at twice the velocity in July as compared with September 
to October (3.6 versus 1.7 mi/day) . Hence, the outer Bay 
regions are colder while the inner or nearshore regions 
continue to warm from atmospheric heating. 

The mean sea surface temperatures in the north- 
ern and southern areas of the Bay are consistently warmer 
than those in the center area for all months except 
November through January (when generally uniform tempera- 
tures prevail) . The northern area exhibited consistently 
higher temperatures than did the southern area throughout 
this period. The warm spot in the northern region is pro- 
bably due to the shallower depths of water there (see 
Figure 12) thus permitting greater insolation and reflective 
heating of the waters as earlier mentioned. Also, the shal- 
lower northern region is somewhat isolated from the effects 
of upwelling and thus does not undergo mixing with the 
deeper and colder waters. 
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3 . Sununary 



Inspection of the mean isothermal surfaces por- 
trayed in Figures 19-30 reveal several blank spots in the 
data field, especially in the northeastern section of the 
Bay. It would be of considerable benefit to conduct a 
further search for data within these blank spots (Blocks 
16, 18, and the non-numbered block shown crossed out in 
Figure 1) thus providing a more complete picture of the 
sea surface temperature distribution. 
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IV. USE OF ISOTHERMAL SURFACES AS INDICATORS OF MEAN CURRENT 



If mean flow can be approximated in Monterey Bay using 
mean isothermal surfaces, then the topographies of these 
surfaces must be associated with all dynamic forces acting 
to establish mean flow. Some of the factors involved are 
as follows: 

A. POSSIBLE GENERATING FORCES 
1 . Offshore Currents 

The offshore currents are considered by some to be 
the major driving forces behind the generation of Bay cur- 
rents. Some theorists believe that a northerly flow adja- 
cent to the Bay will establish a clockwise gyre within it, 
and conversely for a southerly flow. This theory was pre- 
sented recently by Garcia [1971] using a mathematical model 
of the Bay. 

a. The California Current 

Transgressing a large clockwise circle in the 
Northern Hemisphere, the California Current makes its way 
south along the coast year-round. The southern flow of its 
outer boundaries is felt at the seaward edge of the Bay from 
mid-February until early November, with mean velocities 
ranging from 1.0 to 6.7 mi/day (see Chapter III, paragraph 
E.l). 

b. The Davidson Current 

Prevailing from November to early February, the 
Davidson Current flows north adjacent to the coast, bringing 
relatively warmer waters into the Bay (see Figure 31) . 
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2. Winds 



Long acting and relatively steady winds will 
generate related current patterns wherever and whenever 
such wind systems exist, 

3 . T ides 

Tidal currents are active in Monterey Bay but are 
not considered in this study since they are very nearly 
(i£ not entirely) averaged out in taking monthly means. 

4 . Upwelling 

Upwelling has been studied as thoroughly as any 
oceanographic feature, but almost all studies considered 
only its vertical aspects. As discussed in Chapter III, 
paragraph B, the author found that the horizontal dif- 
ferences in upwelling in the vicinity of the Submarine Can- 
yon were apparently a contributing factor in the resulting 
shapes of the mean isothermal surfaces (possibly a major 
factor during the months of maximum, rates of up/dovmwelling) . 
The possible influence of up/downwelling upon mean indicated 
currents is discussed in paragraph D of this chapter. 

B. VERIFYING THE CORRELATION BETWEEN ISOTHERMAL AND ISO- 

PYCNAL SURFACES 

To employ an approach similar to that used by Leipper 
[1970] in the Gulf of Mexico (discussed in the Introduction) , 
it was first necessary to show that changes in mean density 
were directly proportional to changes in mean temperature 
throughout the water column in the Bay. There was not 
enough data available on mean density to allow the con- 
struction of vertical density sections, therefore, annual 



37 



cycles of mean density and temperature (for several blocks) 
were compared. The results showed that mean density and 
temperature were closely related throughout the Bay and at 
all depths considered except in the region above 10m where 
the correlation was not quite as good. The below listed 
comparisons demonstrated the results: 

(1) Figure 32 (mean density structure, Blocks 1 and 5) to 
Figure 14 (mean temperature structure, Blocks 1 and 5), 

(2) Figure 32 (mean density structure. Block 20) to 
Figure 17 (mean temperature structure. Block 20), 

(3) Figure 33 (mean temperature structure with mean density 
structure (dashed lines) superimposed. Block 13). 

The relatively small total change in mean density through- 
out the v/ater column was not unexpected. Previously measured 
currents in the Bay were generally on the order of U.2S to 
0.5 kt , and rarely (if ever) exceeded 1 kt. Therefore, the 
small changes in density (assuming isopycnal surfaces are 
related to indicated mean current) would imply small values 
of mean current. This appears to be the case. 

C. PORTRAYAL OF THE ISOTHERMAL SURFACES 

The foregoing results provide evidence demonstrating 
the close correlation between mean isothermal and isopycnal 
surfaces. Assuming that Leipper's approximate geostrophic 
approach could be extended from deep waters to the waters of 
Monterey Bay, the mean isothermal surfaces would then ef- 
fectively replace the mean isopycnal surfaces and indicated 
mean currents could be deduced from them. On this basis. 
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topographies of the mean monthly isothermal surfaces were 
assumed to represent mean flow and were portrayed in Fig- 
ures 34-48. 

1 . Selecting the Isothermal Surface for Portrayal 

The annual variation in upwelling, causing the 

isotherms to rise and fall through a considerable depth in- 
terval, prevents the use of any one isothermal surface for 
the whole year. Instead, isotherms were chosen for a given 
month by inspection of selected vertical temperature sec- 
tions. The isotherm that best described the temperature 
characteristics of the mid-waters (40-100m) was chosen from 
the corresponding vertical section. When a shift was made 
from one isothermal surface to another, both surfaces were 
drawn to assist in observing continuity. The mid-waters 
were chosen for two reasons: (1) rapidly changing thermal 

conditions in the surface waters could be misleading since 
they did not correlate as well to mean density conditions 
as did the deeper waters, and (2) the 40-100m region should 
respond to offshore currents, upwelling, and wind (to a 
lesser extent) and thus provide mean isothermal surfaces 
most likely to be related to changes in mean current. 

2 . Description of Figures 54-48 

The comments in paragraph III. E.l concerning loca- 
tion of data points, mean winds, and mean offshore currents 
are applicable to these figures as well. The number of 
temperature values used to compute the mean temperatures 
for each block are listed in Table I (the values at 50m are 
representative of those throughout the mid-waters) . 
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Topographies of the mean isothermal surfaces were 
contoured at 10m intervals. Dashed lines indicate probable 
extensions of contours where data was not available. Depth 
contours of isothermal surfaces are portrayed only in regions 
of the Bay where depths of water are greater than or equal 
to the values of the contours. The termination of a given 
contour is shoA^^n by a short segment of the equal valued bot- 
tom contour drawn across it. 

The letter H (high) indicates the greatest depths 
of the isothermal surface and the letter L (low) indicates 
the shallo\\rest depths. 

D. MEAN CURRENTS INDICATED BY THE ISOTHERMAL SURFACES 

In considering the indicated mean currents it should be 
emphasized that the results are based on the combination of 
data compiled from 1929 through 1968. The various regions 
of "highs” and "lows" appearing on the mean monthly iso- 
thermal surfaces (Figures 34-48) cannot be considered as 
fixed. Skogsberg [1936] showed that internal waves and ad- 
vection, acting in the southern portion of the Bay, caused 
daily (and even hourly) changes in the thermal structures. 

For any given month the contours vary from day to day and 
also from year to year, and therefore, the figures are use- 
ful as indicators of the average thermal conditions and not 
the conditions for any given time. 

Since it was assumed that the topographies of the mean 
isothermal surfaces represented mean flow (see paragraph C) , 
the following rules apply: (1) mean circulation is counter- 

clockwise about regions marked with an L (cold water, high 
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mean density regions) , and clockwise about regions marked 
with an H (warm water, low mean density regions), (2) mean 
flow follows the contours with the colder water always to 
the left, and (3) stronger currents are indicated by closer 
spacing of the contours, and conversely. 

The indicated mean currents (inferred from Figures 34- 
48) are discussed below. To prevent repetition, the term 
"mean flow” will be used to imply indicated or inferred mean 
flow and the term "intensity” (or "strength”) will imply the 
relative intensity based on the spacing of contours. 

1 . Intensities of Mean Flow in the Bay and Offshore 
The intensity of mean flow (I) in the Bay is rep- 
resented by a formula relating it to contour spacing as 
follows : 



I = AD/5L 



where , 

I = Intensity of mean flow in mi/day, 

AD = The difference in depth (m) between contours, 

L = The measured length (in) between contours. 

The factor of 5 (in the denominator) was applied to make the 
annual range of I agree with the annual range of mean off- 
shore current (mi/day) as found in H.O. 507 (see Figure 31). 
The values of I plotted in Figure 31 represent the maximum 
intensities (minimum spacing of contours) for a given month. 
Thus, the assumption is made that the strongest monthly aver- 
age currents in the Bay are equal to the average currents off- 
shore. Months where two isothermal surfaces are portrayed 
are represented by an average value of I . 
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The annual cycle o£ the mean offshore current 
(dashed line in Figure 31) was divided into two periods. 

The California Current prevailed from February (mid-points 
of the months are plotted on the abscissa) until early 
November, the Davidson Current prevailing the remainder of 
the year. The annual cycle of mean current intensity in 
the Bay shows peaks in December, March, and August. The 
peak in December and the general shape of the curve from 
November to February is similar to that of the offshore 
Davidson Current. The peaks in March and August do not ap- 
pear to be related to peaks of the offshore current; in 
fact, the shapes of the Bay and offshore curves approach 
that of a mirror image during the period from March until 
.August. These pe^^ks appear to be related to upwelling and 
are discussed in subparagraph 2 below. 

2 . Annual Variations of Mean F^low in the Bay 

Inferred mean flo^^^ is depicted by arrows with the 
intensity of flow (I) indicated to the nearest whole number 
(mi/day) in Figures 34-48). 

a. November to February (Figures 47, 48, and 34) 

The variations in intensity of flow (as noted 
above) are similar for the Bay and the offshore Davidson 
Current prevailing during this period. As described in para- 
graph A.l, a north- flowing offshore current is considered by 
some to generate a clockwise driven current within the Bay. 

In this case, the value of I (shown in the Figures) 
represent relative intensities at given locations and not 
the maximum monthly values . 
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Referring to the mean isothermal surfaces, the current 
patterns exhibit this clockwise tendency throughout the 
period, suggesting that the Davidson Current might be the 
primary generating force. 

In November, the mean flow enters the Bay from 
the south and spreads into the northeast region with a broad 
clockwise movement indicated in the center region. Flow 
then exits from the northern region of the Bay (note the 
40m contour) . 

In December and January, the same general 
clockwise tendency of flow is indicated. However, the flow 
appears to enter from the northern region and exit from the 
southern region of the Bay. The difference in entry points 
could be due to the offshore current intensity which in- 
creases from 1 mi/day in November to 3.4 mi/day in December 
and January. 

b. February to April (Figures 35-37) 

In February, the "high” (located near the cen- 
ter of the Bay) intensifies and the resulting intensity of 
flow doubles over that in January. Flow appears to enter 
and exit from the northern region following a well defined 
clockwise path throughout the Bay. The clockwise tendency 
of flo\v could not have resulted from offshore current genera- 
tion since the direction of the offshore current is av;ay 
from and nearly perpendicular to the coast (and also at its 
annual minimum speed of 1 mi/day). The maximum rates of up- 
welling occurring in February and March (as discussed in 
Chapter III, paragraph B) may explain the inferred clockwise 
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tendency of the currents. As previously postulated, the 
spreading or up -surge of the cold waters from the Canyon 
into the Shallows regions cause a relatively warm tongue 
(high) in the Canyon region. The high established by this 
up-surge might tend to overshadow the influence of the weak 
flowing California Current during February, 

The 10°C and 11°C isothermal surfaces were con- 
structed for March. Both have similar patterns but the 
gradients are stronger at the greater depths (10°C surface, 
Figure 37) . This further indicates the influence of up- 
welling since upwelling rates are maximum from 60 -100m. The 
currents indicated by the 11°C surface (Figure 36) enter the 
Bay from the south and follow a snake- like pattern before 
exiting from the northern region. The 10°C surface (Figure 
37) exhibits the same general isothermal structure except 
there is another well formed low observed in the northern 
region and the overall intensity of flow is greater. The 
development of the "lows" in March could be caused by an 
influence of the California Current which had increased to 
5 mi/day (near its annual maximum intensity). The influence, 
by the theory of offshore generation, would be to cause a 
counterclockwise gyre within the Bay which would agree with 
formation of a ’’low.'' This was not indicated to be a major 
feature by the isothermal topographies, 

c. August to October (Figures 42-44) 

This period is discussed next because the mean 
flow in the Bay is believed to be a result of conditio3is 
similar to those described above. August and September (the 
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months of maximum rates of downwelling) showed the same dip- 
ping of isotherms across the Canyon and the same relative 
warm spots in the Canyon region as were observed during 
February and March. The possible relationship was discussed 
in Chapter III, paragraph B. The well established highs in 
August and, for the 10°C surface, in September (the deeper 
surface where downwelling rates are also greater) demon- 
strate the possible influence of downwelling. Current trends 
are clockwise (as in February and March) and thus opposite to 
the currents theoretically predicted due to influence by the 
California Current. On the other hand, the 11°C surface 
(shallower waters) in September exhibits more complex trends 
that could have been partly influenced by the California 
Current. Current patterns are counterclockwise in the mid- 
dle and northern regions and a clockv/ise trend is observed 
in the southern region. Downvvrelling is weaker above 50m 
(in September) and apparently does not exert as great an 
influence on the indicated currents. It was also observed 
that flow indicated by the shallow waters (11°C surface) is 
nearly opposite in direction to that indicated by the deeper 
waters (10°C surface) . Since the whole region of the mid- 
waters (40-100m) is assumed to represent the indicated cur- 
rents, it is concluded that the shallow and the deep surfaces 
(in September) combined in some fashion to represent mean flow. 
The best approximation of the mean flow is thus a combination 
of flows inferred from both surfaces. 



45 



d. October (Figures 45 and 46) 

The current trends indicated by the 11°C iso- 
thermal surface are similar to those of the 11°C surface in 
September (Figure 44) . Flow enters the southern region of 
the Bay, follows an "S-shaped” path, and exits from the 
northern region. The high observed in the southern region 
appears to be a holdover from the high observed during July 
and August . 

The 12°C surface (Figure 46) is generally 
similar to the 11°C surface except the high is absent. The 
general counterclockwise trends (expecially for the 12°C 
surf ace) conform' to the theory of generation by an offshore 
current. 

In general, it appeared that weaker offshore 
currents came into the Bay, while stronger currents caused 
a reverse eddy within the Bay. 

e. April to August (Figures 38-41) 

The intensity of indicated mean flov; is low 
throughout this period. Figure 31 shows the weak mean cur- 
rents which decrease to an annual minimum in June. 

Current trends in April are similar to those 
in March (Figures 36-37) but the intensity of flow is only 
half as great. The weakening of the high in the center 
region corresponds to the decrease in the rate of upwelling 
which ends altogether in May to June. The low in the southern 
region remained the same strength as it was in March but in- 
truded further into the Bay (as noted by the 40m contour in 
Figure 38) . 
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The indicated mean flow during May is almost due 
south, dipping into the Bay and broadly following the shore- 
line. The direction of flow corresponds closely to the di- 
rection of mean offshore current in the central and northern 
regions and tends seaward in the southern region. 

A high develops in June which extends through- 
out the Bay. Mean flow enters the Bay from the north and 
follov\rs a broad clockwise path, exiting also from the north- 
ern region. Intensity of mean flow is at its annual minimum 
in June . 

The high persists in July and a low is seen to 
develop in the northern region. Flow enters the Bay from 
the north and is observed to exit from the southern region. 

The development of the high in June and its 
subsequent strengthening from June through August (Figure 
42) could have been partly due to downwelling which com- 
mences in June and ends in early November (see Chapter III). 

In summary, the inferred mean flow in the Bay 
is based solely on the topographies of mean monthly iso- 
thermal surfaces. Inferred current trends could be related 
to generation by offshore currents only during the Davidson 
Current Period and possibly in the month of May. Trends of 
inferred’ current during the remainder of the year appear to 
be influenced more by up/downwelling than by offshore cur- 
rents. The possibility of up/downwelling exerting a greater 
influence than offshore currents on mean flov/ is difficult 
to test and verify. However, the maximum intensity of mean 
offshore current during any month of upwelling or downwelling 



47 



is only 5 ni/day (about l/5th of a kt) , and the annual maxi- 
mum reaches only 6.7 mi/day. Thus, it would appear that 
the relatively small horizontal density changes (assumed to 
be mainly a result of the surging of the cold waters during 
up/downwelling) betxs'een the waters of the Canyon and those 
of neighboring shallow areas could be equally (if not more) 
responsible for influencing mean Bay currents. The above 
phenomena are of course, related in some complex manner. 

There is insufficient data concerning pre- 
viously measured currents in the Bay to verify the influence- 
of up/downwelling on Bay currents. However, the author 
found several limited current studies which could be com- 
pared to the inferred Bay currents of this report. 

E. COMPARISONS OF INFERRED MEAN CURRENTS TO PREVIOUSLY 

MEASURED CURRENTS 

Comparisons between the inferred mean monthly flow 
(based on many years of data) and several short term (6-48 
hours) current measurements during a given month cannot be 
expected to show more than broad agreement for specific 
times. The positions of the highs and lows (and hence the 
distribution of indicated current) on the mean isothermal 
surfaces for a given month represent the most likely con- 
ditions over a long time period and not those for a given 
day of the month, or even for a given year for that month. 
Finally, the short term currents must be considered only as 
general trends since their directions and speeds are heavily 
influenced by tidal forces and local winds. The following 
comparisons are more encouraging than the above limitations 



suggest. 
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1. March 



Stevenson [1954] deployed three surface drogues in 
early March from a location 1/2 mi offshore in the southern- 
most area of the Bay. All drogues moved shoreward (until 
beaching) as shown by the dotted arrows in Figure 36. The 
shoreward movements of the drogues would appear to agree 
with the inferred currents which tend shoreward in the 
region where the drogues were deployed. 

2 . J uly 

McKay [1970] conducted current measurements in late 
July using a Geomagnetic Electrokine tograph (G.E.K.). The 
measured ciirrents (indicated by dotted arrows in Figure 41) 
did not fully agree with the plotted indicated current pat- 
terns: however, the position of the "low” in the northern 
region would support such current movement. The shape of 
the "low" could not be determined since there was no data 
for adjacent blocks to the north and east. 

3 . August 

McKay conducted current measurements with the G.E.X. 
in early and mid-August, 1969 (solid short arrows in Fig- 
ure 49) and Stoddard [1971] employed surface drogues in mid - 
August, 1970 (dotted arrows). 

a. McKay's G.E.K. Measurements 

The shorter solid arrows represent individual 
G.E.K. measurements. Considering the region bounded by the 
60m and 80m contours, the correlation of the measured cur- 
rents to the inferred mean currents is very close, with 
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about 3/4ths of the G.E.K, measurements observed to follow 
the contours. In the region bounded by the 80m contour, the 
measured currents broadly exhibit a clockwise motion with a 
general confusion in direction noted (which agrees with the 
plotted pof'.ition of the high) . 

The G.E.K. measurements to the south of the 60m 
contour appear somewhat confused in direction when considered 
as a whole. The arrows with tails represent measurements 
taken on 5 August while all other solid arrows represent 
measurements taken from 11-12 August. During 11-12 August 
the measured currents in the extreme southern region are 
seen to move northward and then to turn towards the west. 

This trend would appear to agree with the mean flo\\' which 
"tHc s t^v^s c I'* 50'^ 3.nci COr*! fTon? 

southern region and could force the northerly current to 
turn westward. The currents m.easured on 5 August tend to- 
wards the southwest which could be due to the conditions 
such as displayed by the isothermal surfaces in July (note 
the 30m contour in Figure 41) . Since all surfaces are as- 
sumed to represent the mid-month conditions, the mean flow 
in July might influence that in early August. The above 
discussion is conjectura] since insufficient data is avail- 
able to construct contours in the region of interest. The 
intention is to show that the current trends by McKay in 
this (southeastern) region show plausible agreement with 
those inferred from the isothermal surfaces. The measured 
current velocities were all significantly higher than inferred 
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values (I) . This is as expected since they were virtually 
instantaneous measurements . 

b. Stoddard's Drogue Study 

The longer dotted arrows in Figure 49 represented 
surface drogue tracks. Six drogues vvrere tracked for about 17 
hours and one (the southernmost shorter track) was tracked 
for 7 hours . 

The two tracks in the region bounded by the 70m 

contour follow definite clockwise gyres, not too far removed 

* 

from the location of the high (generally "agreeing with its 
plotted position) The —remaining tracks do Ifi’ot correlate 

with mean flow; in fact the three tracks to the south are 
nearly 180 degrees out. Data was not av'ailable during August 
in Blocks 10, 11, 14, and 21S (see Figure 1) and therefore, 
the isothermal surfaces (as shown for this region) could be 
in error. 

4 . September 

Stoddard tracked six surface drogues in late 
September, 1970. The tracks show good correlation to the 
shallower 11°C isothermal surface (Figure 44) but generally 
are not well related to the 10°C surface (Figure 43) . The 
following discussion is based on comparison to the 11°C sur- 
face (although the best approximation of indicated mean flow 
would lie somewhere in between the two surfaces) . 

The duration of the drogue tracks varied from 12 to 
25 hours. The author computed the speed of each drogue based 
on the distance made good from deployment until the end of 
tracking. The average value of the six speeds ^^^as 2.7 5 mi/day 
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which compared favorably to the arbitrary value of I for 
September of 3.6 mi/day (Figure 31). The value of I rep- 
resents the maximum mean intensity (minimum contour spacing) 
for a given month, and was made to agree (in total range of 
values) with the intensity of mean flow for the offshore 
current (see paragraph D.l). 

Drogues 1-3 (see circled numbers adjacent to dotted 
arrows in Figure 44) all exhibit clockwise trends which 
broadly agree with the location of the high in the south- 
eastern region of the Bay. Closer correlation is not pos- 
sible since no data exists in the block to the right of the 
high (shown crossed out in Figure 1) . 

Tracks of drogues 4 and 5 cut directly across the 
40m contour; however, their tendencies to move to the right 
when seaward of Point Pinos agree with the general mean cur- 
rent patterns (note the north- flowing mean current which 
would force the currents indicated by tracks 4 and 5 to the 
right) . 

The track of drogue 6 agrees with indicated mean 
flov/, following the 40m contour and tending in a slight 
clockwise path. The speed of. drogue 6 (computed as above) 
was computed to be 4.1 mi/day versus an arbitrary value of 
I equal to 4 mi/day in the same region. 

5 . October 

Stoddard tracked four surface drogues in early 
October, 1970. The tracks show excellent correlation to the 
shallower 12°C surface (dotted arrows in Figure 46) and fair 
correlation to the 11°C surface (Figure 45) . Drogue 3 was 
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tracked for 10 hours and all other drogues were tracked for 
only six hours. The average computed speed for drogues 1, 

2, and 4 was eight mi/day, while drogue 3 made good a speed 
of 6 mi/day. The speeds inferred from contour spacing are 
only 2 to 3 mi/day. Since the duration of tracking was 
short, the effects of the tidal currents (and local winds) 
could have introduced considerable influence on the speeds 
and directions of the drogue tracks. 

All drogue tracks broadly followed the 30m contour 
(Figure 46) and all except drogue 3 moved in a direction 
related to the indicated mean current. The path of drogue 
4 almost exactly agrees with the indicated flow. 

6 . November 

Stoddard tracked eight surface drogues in mid- 
November, 1970 (shown by solid, dashed, and dotted arrows, 
numbered 1-8 in Figure 50). The duration of tracking varied 
from 6-42 hours. The computed speeds (again based on dis- 



tance made good per 


total tracking time) 


are listed below: 


Drogue No. 


Duration of Tracking 
(hr) 


Computed Speed 
(mi/day) 


1 


.42 


5.7 


2 


22 


3.9 


3 


21 


1.2 


4 


23 


3.4 


5 


7.5 


11.2 


6 


6.5 


14.1 


7 


6.5 


10.7 


8 


6 


10.8 



The consistently greater speeds for drogues 5-8 
were probably due to tidal currents- prevalent in the short 
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duration of tracking time and were thus disregarded. Drogues 
1-4 exhibited speeds closely related to the arbitrary value 
of I for November. The combined average speed of drogues 
1-4 was 3.6 rai/day while the arbitrary maximum value of I 
was 4.2 (Figure 31). The average value of I for November 
would be somewhat lower (but nearly impossible to estimate 
from contour spacing alone) . 

The tracks of all drogues exhibited clockwise 
tendencies which correlated with indicated mean current 
about the extensive high in the center of the Bay. The in- 
dividual tracks of drogues 1 and 8 closely approximate the 
indicated mean current while the remaining drogues showed 
fair correlation (Figure 50) . 

It is noted that the speed of drogue 3 (1.2 mi/day) 
is significantly lower than the speeds of drogues 1, 2, and 
4. Also, the track of drogue 3 describes a tight clockwise 
path which is not the case for the other tracks. This led 
the author to believe that drogue 3 might have been nearer 
the center of a high than is shown in Figure 50. It is as- 
sumed that the high (and the entire mean isothermal surface 
as well) could probably vary in location from day to day 
(and year to year) for a given month. On this basis, the 
mean isothermal surface for November was translated to a 
position where the high overlayed the track of drogue 3. 

The distance of translation was 3 mi in a west- southwest 
direction. Figure 51 shows the results of the translation. 
All drogue tracks are now observed to closely follow the 
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translated contours and the correlation to indicated mean 
current is excellent. 

F . SUMMARY 

The general agreement between measured currents from 
past studies and indicated mean currents for this study is 
encouraging. The total number of comparisons made is far 
from that required to fully substantiate the inferred mean 
currents of this study; however, the comparisons do lend 
considerable support to the feasibility of using mean iso- 
thermal surfaces as indicators of mean flow, which was a 
primary objective of this study. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

1. Progressive warming with time is indicated through- 
out Monterey Bay, mainly occurring at depths from the surface 
to 300m. The temperature was noted to increase as much as 
1.8°C over the 29 year interval from 1931 to 1960. 

1 . The intensity of upwelling varies in different geo- 
graphic regions of the Bay. This variation in the maximum 
rates of upwelling (and downwelling) is believed to be a 
major factor controlling the shapes of mean isothermal sur- 
faces (and indicated mean currents as well) during the months 
of maximum up/downwelling . 

Wk /-N-C yrf\ ^ rk 1 ^ • -v* C* 

% XJi^ \J JL. KA. ^ i.L ^ Cl-O 

indicators of mean current shows promise. The comparisons 
of inferred mean flow to direct current measurements (as 
shown in Figures 34-48) demonstrates general agreement be- 
tween the two. 

B. RECOMMENDATIONS 

1. The data collected throughout Monterey Bay is in- 
sufficient for several of the blocks in the author's grid 
system (Figure 1). To obtain a more complete picture, both 
in the mean monthly sea surface temperature distribution 
and in the topographies of mean monthly isothermal surfaces, 
a continued search of data for the following blocks would be 
beneficial: 7, 8, 9, 10, 11, 14, 16, 18, 21S as well as the 

block shown crossed out in Figure 1. 
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2. Future studies should be made of the mean monthly- 
current distribution in the Bay (the author was unable to 
locate any previous work in this area) . Whatever methods 
are used in the future might be comj)ared to the mean flow 
inferred in this study for purposes of attempted verifica- 
tion. 

3. To aid in the future study of oceanographic para- 
meters (in general) of Monterey Bay, a computer storage and 
retrieval system should be instituted and maintained on a 
continuing basis. An example of a computer system designed 
for this purpose is presented in Appendix B (the "final 
computer program, number 11"). 

Closer coordination with NODC is required of all 
oceanographic facilities so that the oceanographic data 
collected in Monterey Bay over the years can be made avail- 
able to any investigator desiring its use. 
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APPENDIX A; DATA SOURCES 



A. PROBLEMS ENCOUNTERED 

The only previous research found that categorized oceano- 
graphic data collected in Monterey Bay was that presented in 
an NPS research paper by P. M. Magruder, Jr., [1969]. 
Magruder's paper was most beneficial in delineating the major 
sources of data. Many reports were located covering related 
subjects such as biology, sound velocity variation, and 
light attenuation studies. Many of these reports used tem- 
perature and salinity data. However, the data was generally 
not contained therein, nor is there any reference to its lo- 
cation in some reports. There is a general lack of continu- 
ity in repoi Ls waicli contained temperature and salinity data. 
Methods of collection were sometimes omitted and very seldom 
was there any information describing environmental condi- 
tions, etc. 



B. MAJOR DATA SOURCES (FIGURE 3) 

1 . Bolin Data 

Rolf L. Bolin and collaborators [1964] collected 
Nansen samples over a five year period from 1951-1955 at the 
position 36-42N, 122-02W. Casts of 24-bottles were gener- 
ally employed to depths of 1000m at weekly intervals. A 
total of 233 casts were made and analyzed for temperature, 
salinity, and other chemical parameters. A copy of this re- 
port is on file at the Naval Postgraduate School Library. 
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2 . FI eet Numerical Weather Center (FN'.VC) XBT Data 



Nansen cast, bathythermograph (BT) , experimental 
BT (XBT), and message BT data are on file at FNWC. The 
Nansen cast and .BT data are also filed at the National 
Oceanographic Data Center (NODC) and are kept more com- 
pletely there. 

FNWC is presently digitizing the XBT data for in- 
clusion in computer tape files at NODC. Two 9 -track com- 
puter tapes (NFS 264 and NFS 265), lield by Dr. Denner at 
the Naval Fostgraduate School Oceanography Department, 
contain the world oceans XBT data which consists of about 
40,000 individual stations. A computer program (see article B.2, 
Appendix B) written by the author reveals only 43 XBT sta- 
tions scattered randomly throughout the Bay. There have 
been considerably more than 43 XBT drops in the Bay over the 
last three years but the data is not as yet available. 

3 . Hopkins-CalCOFI Data 

Hopkins Marine Station has conducted a continuous 
hydrobiological survey of Monterey Bay for the California 
Cooperative Oceanic Fisheries Investigations (CalCOFI) since 
1951. Commencing in March of. 1954 the data (unpublished) 
was tabulated and distributed in annual reports containing 
monthly averages of Bucket, Nansen, and BT temperatures, 
salinities, and densities (sigma-T) for six oceanographic 
stations in the Bay. The averages were based on approxi- 
mately v.'eekly visits to each station; [Hopkins 1954-1967]. 
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a. Ocean Temperatures 

Surface temperatures were read from a bucket 
thermometer and BT data was listed at 10m intervals from 
the surface to 50m for stations 2, S, 4, and 6, and to 30 
and 20m for stations 1 and 5, respectively. This data was 
not computerized since the majority of it was included as 
NODC-BT data in paragraph 4 below. 

b. Salinity and Sigma-T Data 

Nansen samples taken at the surface and 15m 
were titrated for chlorinity and thence salinity and sigma-T 
were determined. Temperature, salinity, and sigma-T were 
tabulated along with their monthly averages for each station. 

4 . National Oceanographic Data Center (NODCj Data 
Magruder obtained complete printed listings of 
Nansen cast and BT data from NODC in January of 1969. The 
area of coverage included Marsden Squares 121 (61 and 62). 

The listings are available in the NPGS Oceanography files. 

For purposes of, (1) rapid computer processing, 
and (2) updating of data through 1970; magnetic tapes were 
requested from NODC. Two 7- track tapes ;>rere received in 
December 1970, one listing Nansen cast data in Marsden 
Square 121(R1019) and the other listing BT data (R1030) . 

The formats of these tapes were different and had to be pro- 
cessed separately (see articles B.4, 5 of Appendix B) . The 
tapes are the property of the NPGS Oceanography Department 
(see Dr. Leipper) . 
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The total number of Nansen and BT stations within 
Monterey Bay (retrieved from above tapes) were 107 and 1377 
respectively . 

a. Nansen Cast Data 

The data included temperature, salinity, oxygen, 
phosphates, silicates, nitrates, and pH as observed para- 
meters. All values were interpolated to standard depths by 
a three point Lagrange method. Each station is identified 
by a Master Card, which lists such information as ship name, 
position, date, etc. 

b. BT Data 

BT temperatures were listed at 5m intervals 
from the surface to maximum depth. Each station was again 
identified by a Master Card. 

5 . Naval Postgraduate School (NPGS) Data 

Four stations were established by NPGS in October 
1966. Nansen casts and BT drops were made regularly and 
all chemical parameters were filed quarterly. 

A continuing program of computer processing and 
filing of data cards for these stations is in progress. 

6 . Skogsberg Data 

Tage Skogsberg [1936] conducted an exhaustive study 
of thermal conditions in the southern portion of Monterey 
Bay, collecting over 15,000 samples from 1929-1933. Data 
included temperature, chlorinity, phosphate, and silicate 
from 23 separate stations. The original data is on file at 
Hopkins Marine Station. Skogsberg tabulated monthly temperature 
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averages for two of his stations, both of which were used in 
the author's analysis (Chapters III and IV). 

7 . Sea Surface Temperatures 

a. Hopkins Marine Station 

Sea surface temperature and salinity has been 
measured daily at the station since 1919. All records are 
filed by Hopkins and monthly averages are published by 
Scripps Institution of Oceanography under the title: 

"Data Report, Surface Water Temperatures at Shore Stations,-' 
United States West Coast and Baja California." 

b. Santa Cruz 

Sea surface temperature has been measured daily 
since 1956 and data is also published by Scripps under the 

r*orr»o oo 

^ ^ ^ ^ ^ • 

C. MINOR DATA SOURCES 

The below listed sources comprise data covering either 
small areas or that limited to only several observations. 

A complete listing w'as impossible. 

1 . "Monterey Bay Bibliography" [Baron. 1971] 

This comprehensive study was recently completed by 
the Moss Landing Marine Laboratories' of the California State 
Colleges. It includes all known scientific literature con- 
cerning Monterey Bay. Limited copies are now being published. 

2 . Yeske and Waer Data 

Yeske and Waer [1968] studied light attenuation at 
two stations in the Bay during the period 26 July to 31 
August, 1968. A total of 39 Nansen casts were made and BT 
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drops accompanied eacli cast. Data was filed with NODC and 
is contained on the computer tapes referenced in para- 
graph C.4 above. 

3 . Message BT Data 

Very little data was found at FNWC (only 21 sta- 
tions for all of Marsden Square 121) and none more recent 
than 1968. 

.4 . OC-3520 Data 

The Chemical Oceanography Course (OC- 3520) at .NPGS 
sponsors .limited duration Bay cruises for students. Various 
oceanographic parameters are collected and final reports 
are required. The data is not currently being filed but 
some of it is available. 

D. DISCUSSION 

It was clear that a complete source listing of the 
chemical parameters in the Bay \^:ould be a monumental task. 
NODC maintains a computer listing of oceanographic data 
throughout the world's oceans but little of the data col- 
lected in Monterey Bay has ever reached there. Only 107 
Nansen casts were recorded by NODC for the last 40 odd years, 
representing a trifle of the data actually collected. 
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APPIiNDIX B: COMPUTER PROCESSING OF DATA 



Two major goals were sought in designing the "final 
computer program": (1) to store temperature, salinity, and 

density data in a form suitable for the computation of mean 
monthly averages at given geographic locations, and (2) to 
design a flexible system for storage and retrieval of various 
oceanographic parameters on a continuing bases. 

The second goal was achieved by designing the "station 
title card" in such a manner that a simple extension could 
be made to the "station data card" to include additional 
selected oceanographic parameters without modification of 
existing formats (see paragraphs A. 2. a ^ b) . The institution 
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oceanographic data from Monterey Bay a very simple task. 



A. STORiXGE OF TEMPERATURE, SALINITf AND DENSITY (SIGMA-T) 
DATA 

1 . Problems Encountered 

About 60% of available data was already in computer 
format, either on magnetic tape or on data cards. Regret- 
tably, no two formats were consistent. Thus, a separate 
computer program had to be written to reduce each data source 
to the author's format. This often included calculation of 
densities (when they were not included as part of the 
original format) and alv-;ays required interpolation of one 
form or another to place the data at the author's selected 
"standard depths" (standard depths are listed in Table IV) . 
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The reduction of data already in computer format was actu- 
ally more time consuming than processing the tabulated data. 

The remainder of the available data was tabulated 
in the various reports cited in Appendix A. First it was 
keypunched onto computer cards in various formats, depending 
on the format of the data source. The next step was the 
printing out of data on computer paper to compare it for 
accuracy against the original source. After proofreading, 
computer programs were v^/ritten to convert the data on punched 
cards to the desired format. Again, this included calcula- 
tions and interpolations for the desired results. 

The last step in the process was to check the 
various outputs for correctness and then write the data onto 
a computer "master tape." Paragraph B describes the various 
computer programs to the extent of including information not 
apparent from inspection of the programs themselves. 

2 . Format of the Master Tape 

The master tape (NPS serial number 306) served as 
the storage location for all data processed. The general 
format of data includes a station title record (identifying 
the oceanographic station) and a station data record (list- 
ing the oceanographic parameters). Table II lists all data 
contained on the master tape. 

a. The Station Title Record 

The station title record allows the flexibility 
required in general retrieval of data. Referring to Table 
III, the "80-colunm record numbers" refer to the inclusive 
positions on the station title record where the information 
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TABLE II. Listing of Specific Data Contained on Computer 
Master Tape, NFS 506 



Data Source 


Data 




Year (incl . ) 




Remarks 


Bolin 


T,S,o^ 




1951-1955 


229 


Nansen stations 


FNWC 


T 




Various 


40 XBT stations 


Hopkins -CalCOF I 


T,S,a^ 




1954-1967 


972 Nansen stations, 
each station repre- 
senting the monthly 
average of 4 stations 
at 0 and 15m. 


NODC 


See Chap 


II 


Various 


107 


Nansen stations 


NODC 


T 




Various 


1377 


BT stations 


NPGS 


See Chap 


II 


Sept 1966- 


175 


Nansen stations 



Sept 1967 



Note 1: See Appendix A for detailed discussion of data 

sources . 

Note 2: In cases where only portions of a data source were 

computerized, the year and type of data defines 
such portions. 

Note 3; The "Data Source" is listed by author or collecting 
activity. 
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TABLE III. The Station Title Record 



80-Column 



Record No. 


Variable 


Definition 


1-2 


NCARDS 


The total number of station data 
records following the title record. 


4-7 


NANS, XBT, 
BT, STD 


A left-adjusted variable defining 
the type of data collected. 


9-10 


NDAY 


The day data collected. 


12-13 


NMO 


The month data collected. 


15-16 


NYEAR 


The last two digits of the year 
data collected. 


18-21 


NTIME 


The time (GMT) data collected. If 
unknown, NTIME appears as 0000. 


24-25 


LATDEG 


The degrees of station latitude. 


27-30 


RLATMI 


The minutes and tenths of minutes 
of station latitude. 


33-35 


LONDEG 


The degrees of station longitude. 


37-40 


RLONMI 


The minutes and tenths of minutes 
of station longitude. 


43-45 


NWF 


NWF=999 specifies Hopkins -CalCOFI 
monthly average values (4 stations 
per average) . 


48-63 




Ship name. 




or, ■ 




48-49 




NODC country code. 


50-51 




NODC ship code. 


52-55 




NODC station I.D. number. 


56-59 




NODC Stano number. 


65-80 


REF 


Reference data by author or facility 
controlling data collection. 
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was entered. ''Variable” was the computer name (applied by 
author) used in programming. Desired data may be retrieved 
from the master tape by any one, or any combination of the 
listed variables. For example, assume an investigator is 
interested in conducting a study of experimental (XBT) versus 
mechanical (BT) bathythermograph data in the Submarine Can- 
yon, collected over the last five years. By a series of 
"computer logical IF statements” [see Blatt, 1968] all sta- 
tion title records are eliminated except those that are 
either XBT or BT , have geographic positions in the Canyon, 
and fall within the time frame (NYEAR) of 1966 or later, 
b. The Station Data Record 

Each station data record contains the standard 
depth of observation (m) plus one or more of the parameters 
listed in Table IV above the dashed line. The oceanographic 
parameters below the dashed line were not included on the 
master tape but are listed as the recommended additions re- 
ferred to in the second paragraph of this Appendix. 

B. COMPUTER PROGRAMS USED IN DATA STORAGE 

Ten computer programs were required to compile the mas- 
ter tape and are briefly described below. The programs will 
be found in the section entitled, "Computer Programs.” 

1 . Bolin Data 

a. Program Number One 

Cards were punched from the tabulated data of 
depth, temperature, and salinity. Program number one cal- 
culated densities and linearly interpolated all data to the 
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TABLE IV. The Station Data Record 



80-Column 



Record No. 


Variable 


Definition 


1-5 


STDEP 


The author's selected standard depths 
(m) , as: 0000. Standard depths used 
were: 0, 10, 15, 20, 30, 50, 75, 100, 

125, 150, 200, 250, 300, 400, ... , 
1400, 1500, 1750, 2000, 2500, 3000, 
4000, .. , 6000. 


7-11 


TEM 


The temperature to hundredths of 
degrees C, as: 10.25. 


13-17 


SAL 


The salinity to hundredths of parts 
per thousand, as: 33.65. 


19-23 


SGT 


Sigma-T (density) to hundredths of 
units , as : 26.77. 


25-27 


XOXG 


Oxygen to tenths of ml per 1. 


29-31 


XP04 


P04-P to tenths of g-at per 1. 


33-35 


XN02 


N02-N to tenths of g-at per 1. 


37-39 


XN03 


N03-N to tenths of g-at per 1. 


41-43 


SI04 


Si04-Si to tenths of g-at per 1. 


45-47 


XPH 


pH of seawater to tenths of units. 


49-80 


— 


Can be used for additional storage 
as required; e.g., computed sound 
velocity, dynamic parameters, etc. 



Note 1; Data above dashed line is included on master tape. 

Data below dashed line is the recommended additions 
for a future storage and retrieval system. 
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selected standard depths. Linear interpolation was considered 
accurate since the niaximum difference between tabulated and 
standard depths was only 25m. The final output consisted of 
computer cards in the format for the master tape. 

2 . Fleet Numerical Weather Center (FNWC) Data 
a. Program Number Two 

The existing 9-track tapes (NPS 264 and 265) 
containing world ocean's XBT data through 1968 v/ere pro- 
cessed. Data in Monterey Bay was retrieved, linearly in- 
terpolated, and punched out on computer cards in the format 
for the master tape. Linear interpolation v/as exact be- 
cause XBT data occurred at every depth where temperatures 
deviated from a linear path. 

3 . Hopkins -CalCOFI Data 

a. Program Number Three 

Computer cards were punched from the tabulated 
annual reports [Hopkins 1954-1967]. Program number three 
processed the computer cards, linearly interpolated for the 
10m values (a selected standard depth) , and punched new com- 
puter cards in the format for the master tape. 

Since each station represented a monthly average 
value (generally the average of four stations) the data had 
to be processed differently in the final program. A weighting 
factor (NWF=999) was added to each station title record and 
served as a branch instruction to weight this data "four to 
one" over other data on the tape. 



70 




J, 



t 

,:q 

I 





I 



4 . National Oceanographic Data Center (NODC) Nansen 

Cast Data 

a. Program Number Four 

This program read data from a 7-track tape 
(R 1019) , computed the number of station data records fol- 
lowing each station title record (NCARDS) , and wrote the 
data (in modified format) onto a 9-track tape (NPS 297) . 

b. Program Number Five 

Data within Monterey Bay was retrieved from the 
9-track tape (NPS 297) and written directly onto the master 
tape in the author's format. One linear interpolation for 
15m was required per station, all other data being at the 
author's standard depths. 

5 . National Oceanographic Data Center (NODC) Bathy- 
t hermograph Data 

a. Program Number Six 

Data was read from a 7 -track tape (R 1030) and 
written directly onto a 9-track tape (NPS 296). Each sta- 
tion title record contained the number of following station 
data records, hence, no further programming was necessary. 

b. Program Number Seven 

Data within Monterey Bay was retrieved from 
the 9-track tape (NPS 296) and written directly onto the 
master tape in the author's format. No interpolation was 
required as data was entered on the tape (R 1030) at 5m 
intervals . 
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6 . Naval Postgraduate School Data 



a. Program Number Eight 

This program was written by Robert Middelberg 
for the Oceanography Department. He was kind enough to 
allow the author to copy and then modify the program for 
usage herein. Data cards prepared by Mr. Middleberg were 
read by the modified program and densities were computed. 

All data was interpolated to the author's standard depths 
using Middleberg 's "three point Lagrange" method. The 
Lagrange interpolation method more closely approximates 
the normal shape of a temperature profile when data is 
widespread (50-150m between observations) as was the case 
here. Finally, data was written directly onto the master 
in tli 0 nutiiior ^ 5 rGrrria.u . 

7 . Transferring Punched Card Data 
a. Program Number Nine 

Program number nine read the punched output 
from programs 1, 2, and 3 above and wrote it onto the mas- 
ter tape. 

8 . Master Tape Record Count 
a. Program Number Ten 

When working with computer tape it is essential 
to keep an accurate record count of data as it is stored; 
otherwise, errors that may be made would be nearly impos- 
sible to locate. Program number ten was used to keep a 
running account of data stored and a double check for cor- 
rect foimiat. 
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c. 



THE FINAL COMPUTER PROGRAM 



Program number 11 computes the mean monthly temperatures 
and densities (at author's standard depths) throughout the 
Bay, using data from the master tape as its input. Table 
V defines the output as it appeared on the master printout. 
Figure 2 sliows the initial field of blocks used for the 
segregation of data. The numbered blocks on the Figure were 
the only ones containing significant data (at least one 
month with three or more temperatures in the block) . This 
program will also be found in the section entitled, "Com- 
puter Programs." 
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TABLE V: Definition of Computer Variables in the Final 

Output of Program Number 11 



Variable 


Definition 


TEMP (AV) 


The average of all temperatures (weighted 
and unweighted) for a given block, month, 
and depth. 


SIGMA-T (AVG) 


As above for density (sigma-T) . 


NR WTD VALUES 


The number of weighted stations processed 


SUM WTD TEMPS 


The sum of all weighted temperature value 


STDEV WTD TEMPS 


The computed standard deviation of all 
weighted temperature values. 



SUM WTD SIGMA-T As above for density. 
STDEV WTD SIGMA-T As above for density. 



NR UNWTD TEMPS 


The number of non-weiglited temperature 
values . 


SUM UNWTD TEMPS 


THo sun. o£ 3.11 hoh" v/ 0 ightod 
values . 


STDEV UNWTD TEMPS 


The computed standard deviation of all 
non-weighted temperature values. 


NR UNWTD SIGMA-T 


As above for density. 



SUM UNWTD SIG.^^A-T As above for density. 



STDEV UW SIGMA-T 


As above for density. 


WTD TEMPS 


A listing of all weighted temperature 
values . 


WTD SIG^L\-T 


A listing of all weighted sigma-T values. 


UNWTD TEMPS 


A listing of all non-weighted temperature 
values . 


UNWTD SIGMA-T 


A listing of all non-weighted sigma-T 
values . 



Note 1: All variables above are for a given block, month, 

and standard depth. 
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NOTE: Only the nu:;ibered blocks contained significai'it data. 



Figure 2. The Initial Coriputex Grid 



76 




I 



I 



■If 



j 





77 



TIME (MONTHS) 



DEPTH 

(m) 



JFMAMJJ ASOND 




Figure 4 . Average annual temnerature (*C) 
cycle based on five years of data (upper lOOra 
from 1929-1933 f below 100 m from 1932 - 1933 ) by 
Skogsberg. Dashed lines note scale change. 
Note: The 9 ° isotherms by Skogsberg (S) , 
Bolin (B) and CalCOFI surveys (C) demonstrate 
the downward migration of isotherms. 
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NOTE: Dashed line denotes scale change. 



Figure 5. .\imual Mean Temperature Cycle by Bolin 
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NOTE: Dashed line denotes scale change. 



Figure 6. Annual Mean Te;..perature Cycle Based on Ca C0?1 
Surveys . 
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Figure 8. Standard Deviation of 'Ti'emperature , 
Averaged in Typical Blocks, as a Function of 
Depth. 
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Figure 9. 



The Annual lean Tenpcraturc Cycle 



at Blcck 1^. 
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Figure 11. Vertical Mean Temperature S()Ctions S-1 and S-2 for the Month o£ 
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Figure 13. 



The Annual Cycles of Mean Sea Surface Ten,, .r, a . ere 
in Very Shallow regions of tiie Bay. 
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Figure 14. The Annual Mean Tciaperature Cycles at Blocks 
1 and 5 . 
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Figure 15. 



T’.ie Annual Cycles of 'lean Sea Surface Tenpcr a -urc 
in Relatively Deeper Regions of the ShaljC "s. 
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NOTE : 



The steep slopes of the 
Block 4 illustrate the 
waters into the region. 



10°C and 11°C isotherns at 
up-slope entry or cola 
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Typical Annual 
Canyon vei'sus 



'^ean Teripevature Cycles 
an 0r5shore Station. 
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NOTE: The intensity of flo’.v in the Bay (1) is based on 

contour spacing and is shcv:n by a solid line. The 
dashed line denotes the mean offshore flo-v. 



Figure 31, 



Intensities of 'lean FIc in the Bay n' 
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NOTE: The dashed lines 

i'i<*ure 32. Tiie Annual 



denote ;aid- values of ,:iean density, 
ean Density Cycles at Selected Slocks. 
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NOTE; 



ihe non-labeleci dashed lines denote nid-values of 
mean density. 



■figure 33. 



Tiie \nnual i'ean Te.iper ature versus ensit^ 
Cycles at Block 11. 
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WRITE (6, 10 5) STOEPC JJ) tSTEM( JJ) tSTSAL(JJ) ,STSGT(JJ} 
1C5 FORMAT ( 1X,F5.0,1X,3(F5.2,1X)) 






•5^ 



(\i ro lAvOr-oocTNO 

m ro rn nronfnro>^- 



rost 

'TnJ- 



■>> 

■Jf 

•J'r 

•«« 

•jf- 

•jf 



-Ji- II 

















rr 






00 


















rr 






o 
























•• DC 
























Q< 


h“ 






















xo 


< 






















< 


X 






















cc 


cc 










m 






•}$ 






-X 


o 










O 














X 










X 












ooo 












HH 












r-HCXcO 


h-co 










a 


X 










o< 1 


cr^ 








o 


CL 


in 




•If 




r-H 


r-^OO 


(njq:: 


1 — 






o 


CD 


r» 




■If 




r- 


OC h- 


i/.O 


o 






in 


H- 


vO 




•If 




o 


II LU 


CLX 


00 






• 


00 


HH 




■If 




r-H 


(X-JOO 


XI- 








1— 




H 










LL’KX 


II o 


00 






o 


CO 


*• 




-If 




X 


COhh>:_ 


q:< 








• 


H" 


II 




•If 


- 


< 


--h-o 


LU 








00 


LU 


^ < 




If 


00 


"O 


LU ^ 


cox 








1- 


CO 


vO H 




4f 


OC 




Q-XO 


^ HH 






< 


UJ 




HH 00 




^f 


LU 


CvJ 


<oo 


LU 






h- 


00 


< 


•• 




-If 


Z'Z 


(\1 


h-< 


ClO 


-J 




< 




K 


- ^ 




-If 


X 




LUX 


<x 


< 




o 


< 


< 


II < 






< 




iX hh 


h-HH 


(/) 






H 


Q 


HH LU 




-If 


-J 




OO 


00 


h- 




00 


< 




o 






• 




<o*^ 


Xoo 


00 




«• 


O 


O 


•• o 




•If 


o 




ecu. 


OLU 






X 




O 


o 




^f 


• 


00 


H X 


<o 






H- 






X -J 




^f 




cc 


1 OOi — 


CCO 






X 


< 




HH < 




•If 


■w 


LU 


h-QHH 


h-a: 


*~5 




o 


H- 


X 


X H- 




^f 




y~ 


cC-X 


1 a. 


W 






O 


< 


CL o-' 




•If 




X 


x< 


O' 








K 


X 


O ^-co 




■If 


o 


< 


OOQ 


o: 


LU*-. 


w 


q: 


W 


1 — 


K HH 


rH r-H 


•If 


o 


-J 


iX cc 


00 lU 






LU 


OJ 




to •• 


H H 


■If 


LU 


• 


U.<<t00CDh- 


OOX 


(X 


X. 


cc 


LU 


♦. <•. 


nn 


^f 


o 


o 


h-OOCL< 




z: 


1- 


ou 


CX 


OC 1- H 


W-. w* 


■If 




• 


<< 


ct^XO 


. •«. 




o 


X 


o 


O <CO 


H H 


-If 


K 


o 


h-OLU< 


-OC\J 




2: 






cC OQ 




•If 


LL O 




< OOOeX 


-> • 


u 


< 


LU 




cc cc 


>>- 


•If 


O II 




QLLH- 


1 - o 


wLT^ 


HH 




> 


• • 


LU »"<LUO 


OO 


-If 


vfO 




Ohh<2LL 


CLU. 


h- 


H 




a: 


hhOO 


W'W 


■If 


<ro 




X h-HO 


LUw 


< 


UJ 


QL 


a 


O hh LU 


II II 


^f 


o • 




LUjX 




on 




CO 


ct: 


cc 


< •'Qq: 


LULU 


•If 


hLU 


•S^ 


coxrc 


OXIU 




CO 




< 


cc 


LU hh*^ 


OO 


-If 


h-X 




X UJ 


LUC- 


COX 




< 




UJ 


*^CC »-<UJLU 


<< 




CX)HH 




<00X“0h-< 


f-4 


y- 


H- 


a: 


•• 


«• » 


XD. 




vOH 


X 


■ xlu 


l-h- 


«>* 


X 


< 


LU. 


»• 




COCO 


■if 


O - 


< 


KCOCOhh 


nOO 


UJ 


O 


>r-x 


cox O^XCO 




■If 


wO^f 


a: 


oiDujuucx::^: 


o • 


X 




ajCf-4 


O-H 0»-^ 


OO 


■If 


o 


O 


ox h" h* Xo 


r-4CT\ 






r— H X, 


r-HV. »-4X.O 


oo-;t- 


■If 


CCOO 


O 


OX hh 


HH < 


•‘CLOJHUJXUJ 




HH CO 




•If 


Of”0 


fX 


xoo:cl:<ci: 




OLUXQnO*-- 


vO^ -^vO^ 


hH r-H Q 


■If 


-)cx: 


cc 


ox 


XHH- 




-J' 


1 — 


h-'H HH^w| J 


OOQ 


•If 


ox 




CO 


< 1 


UU<»-< 




OLL'<l 


LU< C\JLU<< 


OO 


■If 


LLhh 




o • • 


• OO' 


k5:k 




l-h-t-Ot-Sia-l-SlO- 1 vO II 1-511-0. 


LULLX 


■«■ 


>t CO 




<*-^(\)n-if 




hhcxo»-< II »— »HHCt:oocDor^»--< 


■r- 


o > 




LU 


•If <c 


o:oo: 2 :oo:^Ocoa:oi— or.oi— ii *-'^q^ok-k 2 :oo.’oo 


•if 


OOCO 




OC 


■If ^ 


3:ii-0 


o a 




jCLUCO:SLL00hhhh»-h 3:LL- tOLUJ— 


-If 


OLU e 




















LLU.CO 


■if 


XXH“ 








nOO 




o o 




h- 


coo O O' 


t • • 


•If 


xajoc 








OC\J 


n o 


o 


oo o o 


ooo 


-If 


< o 












CNJ 


r-H 


sThh in hh 


ooo 


•If 


O LU 









o o o 




o 



133 



DIMENSION INnATA(19) , N SA VE ( 3 5 , 1 9 ) 

DATA NT0T/0/,H/0/, JTOT/0/, 1/0/ 

NT0T=T0TAL NR OF 80 COLUMN RECOFOS WR I TT EN ( NT0T= M+ JTOT ) 
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